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Abstract 

Ferromagnetic-insulator (FI) based Josephson junctions are promising candidates for a coherent superconducting 
quantum bit as well as a classical superconducting logic circuit. Recently the appearance of an intriguing atomic-scale 
0-7T transition has been theoretically predicted. In order to uncover the mechanism of this phenomena, we numerically 
calculate the spectrum of Andreev bound states in a FI barrier by diagonalizing the Bogoliubov-de Gennes equation. 
We show that Andreev spectrum drastically depends on the parity of the Fl-layer number L and accordingly the 7r 
(0) state is always more stable than the (n) state if L is odd (even). 
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1. Introduction 

The peculiarity of the proximity effect in 
superconductor/ferromagnetic- metal (S/FM) bi- 
layers is the damped oscillation of the pair ampli- 
tude inside a FM [1,2]. This anomalous proximity 
effect leads to the n Josephson S/FM/S junc- 
tion [3,4] which has the opposite sign to the super- 
conducting order parameter in two S electrodes in 
the ground state. Experimentally 7r-junction was 
firstly observed by Ryazanov [5] and Kontos [6] and 
since then a lot of progress has been made in the 
physics of 7r-junctions and now they are proving to 
be promising elements of superconducting classical 
and quantum circuits [7-10]. 

On the other hands, recently a possibility of 7r 



junction formation in a Josephson junction with a 
ferromagnetic insulator (FI) has been theoreti- 
cally predicted [11-20]. The it junction using such 
an insulating barrier is very promising for future 
qubit [21-24] and microwave [25] applications be- 
cause of the low decoherence nature [26,27]. More 
importantly, it has been shown that the ground state 
of S/FI/S junction alternates between 0- and ir- 
states when thickness of FI is increasing by a sin- 
gle atomic layer [16,18]. In this paper in order to 
understand the physical mechanism of the anoma- 
lous atomic scale 0-7T transition, we will calculate the 
spectrum of the Andreev bound states in such sys- 
tems. Based on this calculation, we will show that 
Andreev spectrum drastically depends on the par- 
ity of the FI layer number L and thence the ir (0) 
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state is always more stable than the (ir) state if L 
is odd (even). 

In this paper we focused on the one dimensional s- 
wave junction with a FI barrier (Fig. 1(a)). It should 
be noted that the qualitatively same result can be 
obtained for two- or three-dimensional cases. 

2. Model 

Let us consider a one-dimensional tight-binding 
lattice of a superconductor/ferromagnetic-insulator 
/superconductor (S/FI/S) Josephson junction with 
L being the thickness or the numbers of the FI lat- 
tice sites as shown in Fig. 1(a). The lattice constant 
is set to be unity. Electronic states in a s-wave su- 
perconductor are described by the mean-field BCS 
Hamiltonian, 

Hs = ~t C n<x C «'<r + ( 2t "Ms) C n* C ™ 

+# E( c it c l+4 c ii+ h - c -)- (1) 

nGS 

Here (y na (c na .) is the creation (annihilation) opera- 
tor of an electron at a site nGS with spin a = ( j" 
or 4- ) and /i s is the chemical potential. The hopping 
integral t is considered among nearest neighbor sites 
and A is the amplitude of s-wave pair potential. 

The energy dispersion in the Bogoliubov-de 
Gennes picture and the spin resolved density of 
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Fig. 1. (a) The Josephson junction with a ferromagnetic-in- 
sulator barrier on the one-dimensional tight-binding lattice 
and (b) transport through Andreev bound states. An elec- 
tron (black circle) is reflected as a hole (white circle) at the 
interface with the left superconductor and the hole is re- 
flected back as an electron on the right interface. The net 
result is a Cooper pair transfer through the junction. 



states (DOS) for typical FIs are shown schemati- 
cally in Fig. 2. Experimental studies as well as a 
first principle calculations indicate that the band 
structure of an oxide ferromagnet La4Ba2Cu2 0io 
(La422) [28-30] and K 2 CuF 4 [31-33] can be de- 
scribed by Fig. 2(b) in which the up- and down-spin 
bands arc located below and above the Fermi energy 
respectively. The exchange splitting V ex of La422 is 
numerically estimated to be 0.34 eV. Since the ex- 
change splitting is large and the bands are originally 
half-filled, La422 becomes FI with a Curie temper- 
ature of 5 K [28]. Another possible candidates for 
the FI barrier are spinels [34,35], e.g., NiFe2 04, 
rare-earth monopnictides [36-39], e.g., GdN, and 
Yttrium iron garnet (Y 3 FesOi2) [40,41]. 

The Hamiltonian of a ferromagnetic layer can be 
described by a single-band tight-binding model [20] 
as 



n,n'£F,<7 



n£F 



n£F 



V 

2 



cx \ f 



where 



Vox = 4t - 



(2) 



(3) 



is the exchange splitting (g is the gap between up 
and down spin bands) and \l is the chemical poten- 
tial (see Fig. 2(a)). If V cyi > At, this Hamiltonian 
describes FI as shown in Fig. 2. 
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Fig. 2. (a)The band structure of a ferromagnetic-insulator 
in the Bogoliubov-de Gennes picture and (b) corresponding 
spin-resolved density of states (DOS). The energy dispersion 
for a hole with spin a is obtained as a mirror image of that 
for an electron with spin a with respect to Fermi energy. 
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3. Andreev bound states and Josephson 
current 

The Hamiltonian can be diagonalized by the Bo- 
goliubov transformation. Due to the the Andreev re- 
flection at S/FI interfaces, the Andreev bound state 
is formed in the FI barrier (see Fig. 1(b)). Wave 
functions of the Andreev bound state decay far from 
the S/FI interface. In what follows, we focus on the 
subspace for spin-t electron and spin- J, hole. In su- 
perconductors, the wave function of a bound state 
is given by 



*i,(n) =$l 




ik* n 



(4) 



(5) 



Here A, B, C and D are amplitudes of the wave func- 
tion for an outgoing quasiparticle, <j> v is the phase of 
a superconductor, 



^=diag(e^/ 2 ,e-^/ 2 ) 



(6) 



with v = L (R) indicates an left (right) supercon- 
ductor, and 



(7) 
(8) 




with 



ft = VE 2 - A 2 . 



(9) 



The energy E is measured from the Fermi energy 
and 



k = — + i cosh 



A 2 - E 2 



1 
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(10) 



is the complex wave number. In a FI, the wave func- 
tion is given by 



*Fi(n) = 



with 
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and /i , ji , gi and g 2 are amplitudes of wave function 
in a FI. 

By applying the boundary conditions, 



*i(0) = *Fl(0), 
*i(l) = *Fl(l), 

* fl (L) = * F i(L) ) 
y R (L + l) = y Fl (L + l) 



(14) 
(15) 
(16) 
(17) 



we can obtain a secular equation for amplitudes A, 
B, C and D. From this equation, we can numerically 
calculate the Andreev levels £j as a function of the 
phase difference (ft — 4>l — 4>r, where j = 1, • • • , 4. 

The Josephson current can be calculated from the 
Beenakker formula [42], i.e., 



(18) 



where / (e) is the Fermi-Dirac distribution function. 
In the case of a high barrier limit, the Josephson 
current phase relation is described by 

Ij(4>)=Ic sin (p- (19) 
Thus we define the Josephson critical current 1q as 

(20) 



Ic = Ij[- 2 



If Ic is negative (positive), then the tt (0) junction 
is realized. 

4. Numerical results 

In this section, we show numerical results for the 
spectrum of Andreev bound states for a conventional 
S/I/S junction and an S/FI/S junction. In the cal- 
culation, we set fi — ^t s ^ 2t, and A = O.Olt. 

Let us firstly consider Andreev bound states in an 
S/I/S junction. Fig. 3 shows the Andreev spectrum 
as a function of the thickness of the insulating bar- 
rier L. Due to the spin degeneracy, we have 2 An- 
dreev levels for a given <f> and L. It is evident that the 
energy minimum is at <f> = irrespective of the value 
of L. So the overall feature of Andreev levels does 
not depend on L. On the other hand, Fig. 4 shows 
the L dependence of the Andreev spectrum for an 
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Fig. 3. The Andreev spectrum in an S/I/S Josephson junc- 
tion for the case of (a) odd and (b) even insulating-layer 
number L. In the calculation, we set g = 0.25t, fj, = fi s = 2t, 
and A = 0.01*. 



S/FI/S junction. The results indicate that the over- 
all feature of the spectrum strongly depends on the 
parity of L and show that the energy minimum of £\ 
for odd L is at (f> = 0, whereas for even L at <\> = n. 

In Fig. 5, we show the Josephson critical current 
for an S/FI/S junction as a function of L. Tempera- 
ture T is set to be 0.01T C <C T c , where T c is the tran- 
sition temperature of a superconductor. The 7r(0)- 
state is always more stable than the 0(7r)-state when 
the thickness of FI is an odd(even) integer. Based 
on the Andreev spectrum (Fig. 4) the reason can be 
explained as follows. At low temperatures, only the 
Andreev levels below the Fermi energy i.e., e\ and 
£2, contribute to Ic- In the odd (even) L cases, the 
7r- (0-) junction is stable because of 
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Fig. 4. The Andreev spectrum in an S/FI/S Josephson junc- 
tion for the case of (a) odd and (b) even Fl-layer number L. 
In the calculation, we set g = t, /i = fi s = 2t, and A = 0.014. 
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Above analysis provides an new physical interpreta- 
tion of the atomic scale 0-tt transition from the view 
point of the Andreev spectrum. 

5. Summary 

To summarize, we have theoretically studied the 
Andreev levels and the Josephson current in S/FI/S 
junctions by solving the Bogolubov-de Gennes equa- 
tion in order to understand the physical mechanism 
of the atomic scale 0-tt transition. A characteris- 
tic and important feature for such systems is that 
the Andreev spectrum strongly depends on the par- 
ity of the thickness of the FI layer L. As a result, 
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Fig. 5. Josephson critical current Iq at T = 0.01T C <C T c 
as a function of the FI layer thickness L for g = t. The red 
(blue) circles indicate the tt(0) junction. 

the junctions show the atomic scale 0-tt transition. 
Our finding suggests a way of understanding the 
physical origin of the atomic scale 0-tt transition in 
ferromagnetic-insulator based Josephson junctions. 
In this paper, we have only considered the Joseph- 
son transport in the low temperature regime, i.e., 
T <C T c . The calculation of Ic in the finite temper- 
ature region and the analysis based on the Andreev 
spectrum (Fig. 4) are important future problems. 
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